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Abstract 

Trading large volumes of a financial asset in order driven markets 
requires the use of algorithmic execution dividing the volume in many 
transactions in order to minimize costs due to market impact. A proper 
design of an optimal execution strategy strongly depends on a careful 
modeling of market impact, i.e. how the price reacts to trades. In this 
paper we consider a recently introduced market impact model (Bouchaud 
et al, 2004, [1 ), which has the property of describing both the volume 
and the temporal dependence of price change due to trading. We show 
how this model can be used to describe price impact also in aggregated 
trade time or in real time. We then solve analytically and calibrate with 
real data the optimal execution problem both for risk neutral and for risk 
averse investors and we derive an efficient frontier of optimal execution. 
When we include spread costs the problem must be solved numerically 
and we show that the introduction of such costs regularizes the solution. 
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1 Introduction 

Optimal execution of portfolio transactions is receiving a growing attention 
in the recent years. Due in part to the increasing automation of the trading 
process, optimal trading is becoming also an industrial standard, not only for 
very large investors. The problem is generically stated in the following way. 
An investor wants to trade (buy or sell) a given number of shares and wants 
to minimize cost. Execution cost has many components, such as fees, taxes, 
investment delays, market impact, timing risk, and opportunity cost [2J. These 
components depend on the execution algorithm and in many cases are stochastic 
variables. This means that their value depends on the specific realization of the 
trading process. Moreover stochastic costs are associated to a risk, which is 
related to the uncertainty of the cost. As usual in finance, this fact entails a 
risk-reward (in this case, cost) tradeoff, which can be solved by choosing an 
appropriate level of risk aversion. 
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In modern electronic markets, the problem of optimal execution can be a 
multi time scale problem. It is not unlikely that a portfolio manager wants to 
trade an order whose size requires an execution over multiple days. In this case 
we can imagine a three scale decomposition of the execution. First, the port- 
folio manager decides how to split the order across the different days, possibly 
dividing it also among different brokers to minimize information leakage. In this 
case there are specific risks, such as those due to overnight returns. Then, for 
each day the trader divides the day in "macroscopic" intervals, for example 5 
or 15 minutes, and decides how much to trade in each of these intervals. Lastly, 
one has to decide how to trade in each interval, specifying the type of orders 
used (e.g. limit versus market orders) and the strategy to follow (for example, 
when to cross the spread if the price moves in an adverse direction). Although 
this division is very schematic and the optimization is not necessarily segmented 
rigidly in this way, this scheme gives an idea of different types of optimization 
that can be used as building blocks to simplify the whole optimization problem. 

In this paper, we will focus on the second level of trading optimization, 
namely we will assume that a investor wants to trade a position throughout 
the day and that she considers the day divided in intervals of several minutes, 
searching the optimal way of splitting the order among the intervals. Moreover 
we will focus our attention to price impact cost only. However, since we will 
assume that the investor trades by using market orders, bid ask spread costs 
cannot be neglected, and in fact we will see that they play a very important 
role in the solution we obtain. 

There is a vast literature on optimal execution of portfolio transactions. The 
seminal paper of Bertsimas and Lo [3] considered the optimal execution problem 
by modeling a linear and permanent impact and by minimizing the expected 
cost. Not surprisingly they found that the optimal solution is to "slice and dice" 
the order, i.e. to divide it in equal parts and trade it at constant rate. More 
recently, Almgren and Chriss [4] introduced two important innovations. First 
beside the (linear) permanent impact, they considered a temporary component 
of the impact, which takes into account the fact that trading with market orders 
implies a liquidity cost. The second innovation was to include risk aversion in 
the optimization. Since traders can be risk averse, they will minimize expected 
cost for a given level of risk (or, equivalently, minimize risk for a given level of 
expected cost). Therefore Almgren and Chris minimized an objective function 
which is a linear combination of expected cost and risk (quantified as the vari- 
ance of cost), and the coefficient in front of the risk term is a measure of risk 
aversion (see Section 2.1). 

A key aspect in modeling optimal execution is market impact. Roughly 
speaking, market impact is the relation between signed trade volume and price 
change [5]. In general this relation is complicated, because it is a function of 
the type of order, the sign of the trade, the volume, the available liquidity, 
etc.. Several empirical and theoretical studies have considered the properties 
of market impact [HI El HI IS] - The two approaches to optimal execution men- 
tioned above [3l [4], and many others published in the literature, assume (i) 
that contemporaneous impact is linear in volume and (ii) that impact has no 
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transient component, i.e. the trades made at a given time have no effect on the 
returns caused by other trades made at a subsequent time. Empirical research 
indicates that these two assumptions are not consistent with real data. The 
impact of individual transactions is a strongly concave function of the volume 
[7J H] . More importantly, it has been recently shown [TJ [3] that market impact 
has a significant transient component. We should stress that the transient com- 
ponent of impact is strongly related to market efficiency In fact order flow is 
highly persistent in time [iff and, if the impact were fixed and permanent, also 
price movements would be strongly autocorrelated in time. We note that it 
has been shown empirically that the persistence of the order flow is mainly 
determined by order splitting, i.e. optimal execution. 

In this paper we will consider the problem of calibrating the optimal execu- 
tion of financial transactions in the presence of transient market impact. From 
a theoretical point of view this problem has been sketched in [5] and considered 
recently in |12j , where it has been discussed in the framework of the impossibil- 
ity of price manipulation. Some of these theoretical solutions will be re-derived 
here. Moreover our optimal solution will be compared with the one proposed 
by Obizhaeva and Wang jT3] as a solution of a model that takes into account a 
specific mechanism of spread resilience. Then we will consider the more realistic 
case where one considers the bid ask spread cost and we will show how the in- 
clusion of spread cost regularizes the solution. Apart from theoretical analysis, 
the paper will be focused on the calibration on real data. To this end we will 
show how to define and calibrate a transient impact model when we consider a 
macroscopic (i.e. non tick by tick) time measure. In this work we consider real 
time and aggregated trade time. 

The paper is divided as follows. In Section 2 we fix the notation and we 
review the optimal execution problem. In Section 3 we present the dataset we 
used. In Section 4 we review the propagator model in trade time, as introduced 
in Ref. pQ. Section 5 discusses how to extend and calibrate the propagator 
model in real time (the Appendix discusses the calibration in aggregated trade 
time). In Section 6 we discuss the solutions of the optimal execution algorithms 
and in Section 7 we calibrate them to the real data. Finally, in Section 8 we 
draw some conclusions. 

2 Optimal execution problem 

In this section we review the problem of optimal execution and we introduce the 
definitions that we will use in the rest of this work. We define two important 
quantities, namely the trading schedule and the execution cost. 

Consider an investor who wants to trade a block of X shares of a stock. If X 
is positive (negative) the investor wants to buy (sell) . We will consider execution 
programs that make use only of market orders. In this way there is certainty 
that the order will be executed, but the bid ask spread cost might be high. We 
assume that the investor wants to execute her trade in a predetermined time 
horizon T, as, for example, one trading day. We divide T into N intervals of 
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length t = j£, and call t% = kr (k = 0, 1, N — 1) the initial times of each 
interval. The value of N is predetermined by the investor. 

We define a trading schedule as the Af dimensional vector v = (vq, V\, vn-i) t , 
where T denotes the transpose. Each element v k is the signed volume traded 
between time t k and t k+ \. Conventionally, the sign is positive (negative) if the 
agent buys (sells). We clearly have that 

N-l 
3=0 

A trading strategy is now easily defined as a rule for determining the trading 
schedule v k with the information available at time tk-i- Note that at this point 
we are not imposing that the sign of each v k is equal to the sign of X. In fact, 
a buy trading strategy can in principle include sell transactions. 

Each trading schedule is associated with an execution cost. In order to define 
cost, we need to introduce a notation for price. The mid price at the start of 
interval k is denoted by P k and its logarithm is p k = In P k . Note that the average 
price at which a share is traded at time k is not P k . In fact commissions fees 
and, most notably, bid-ask spread raise this price for buy trades or lower it for 
sell trades. We therefore define the series P k which represent the effective prices 
at which shares are actually traded at every step. This is similar to what is done 
in Almgren and Chriss approach (4] (see also below). The precise relationship 
between P k and P k will be clarified later on. 

We define the execution cost C(v) as the difference between the total money 
payed or received during the execution (equal to the sum, over all intervals, of 
the volume traded v k and the corresponding effective price Pk) and the initial 
market value of the position. In formulas it is 



N-l N-l 

£ v k P k - XP = £ 

k=0 k=Q 



C(v) = V v k P k - XP = V v k (P k - P ). (2) 



Execution cost so defined is equivalent to the implementation shortfall intro- 
duced in Ref. [13] . For a buy trade it is the difference between the price payed 
and the price that would have been payed in an infinitely liquid market. We 
note that execution costs needs not to be positive, i.e. one may have nega- 
tive costs (gains) from the execution. For example, if during a buying schedule 
the stock price declines substantially, C(v) could be negative because the total 
money paid is less than the initial market value of the position. 

It is important to stress that, since price is a stochastic process, whose 
properties might be affected by the trading schedule, also execution cost is a 
stochastic process. Different realizations of the price process under the same 
trading schedule give different execution costs. For this reason it is important 
to characterize its expected value and its variance, 

E [C(v)] , Var [C(v)] = E \{C{v) [C(v)]) 2 l . (3) 
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The first quantity measures the expected cost, while the variance quantifies the 
risk associated to the execution. 

2.1 The Almgren Chriss execution scheme 

The problem of optimal execution, as we formulated it, has been the focus of 
many studies. One particularly influential work is the 2001 paper by Almgren 
and Chriss [I] . They assume that the price of the stock at step k is equal to the 
previous price plus a linear market impact term and a random shock: 

P fe = /Vi + 9v k -i + rik-i tj ~ IID(0, a). (4) 

This specification for the price impact is also used in Ref. [3J. 

The first innovations in Almgren and Chriss is that they consider the effective 
price Pfc and they model it as 

Pk = Pk+ pv k + sign(ufe) ■ 5/2, (5) 

where the term sign(iifc) • S/2 is the contribution from the bid ask spread S 
and pvk represents a linear temporary impact. In the terminology of [4], the 
temporary impact accounts for the resilience of the limit orders in the book, 
which relaxes back to the steady state after a trade-induced price movement. 
The equation [2] for the execution costs becomes 



JV-l JV-l N-l N-l 

C(v) = VkP~k-XP = + 6v k ) v j + Y1 (signK)<S72 + pv k ) V k 

k=0 k=0 j=k+l k=0 

(6) 

and the expected value of the costs is 



N-l N-l 

E [C(v)] = -A 2 + + 5/2 ]T \v k \ 



(7) 



fe=0 



fe=0 



If one assumes that all v k have the sign of A, the last term becomes equal to 
XS/2. The optimal solution, i.e. the one that minimizes the expected impact 
costs, is 



'A A A 
v* = argminP fCYvll = ( — , — , — 
6 -- [ N N N 



(8) 



showing that the solution simply consists in trading at a constant rate over the 
periods. 

The variance of the execution cost is 



Var [C(v)] = E (C*(v) - E[C(y)})' 



E 



N-l N-l 



(9) 



We assumed the r]k are independent, so we have E[rjir}j] = for i ^ j and thus: 



JV-l N-l 

Var[C(v)]=^Y / (T, v ^- ( 10 ) 

fc=0 j=k 

The second innovation in Almgren and Chriss is that they consider risk 
aversion for optimal execution: they minimize the sum of expected cost and 
costs' risk. By mimicking the theory on portfolio optimization [15, 16 , Almgren 
and Chriss consider as optimal trading schedule the solution of 

argmin (E[C(v)] + \Var[C(v)}) . (11) 

V 

where A is the coefficient of risk aversion. The higher is the A, the more impor- 
tant is risk with respect to cost. A risk neutral investor corresponds to A = 0. 
The set of solutions to this problem for different values of A is called optimal 
frontier. 

By using the impact model of Eqs. [4] and [5] one obtains 

v k = Acosh(f3(T-t k )), (12) 

where A is a normalization constant. The coefficient f3 is related to the coefficient 
of risk-aversion A, and is given bjQ 



A > 0. (13) 




The solution of Eq. 12 shows that the more risk averse is the investor (i.e. 
the higher is A), the more front loaded is the trading schedule. This means that 
more volume is traded at the beginning of the execution in order to minimize 
the uncertainty on execution price of the last part of the trading schedule. 



3 Data 

For our empirical research we have used two different datasets. The first one 
reports the trading activity at the London Stock Exchange (LSE) and covers 
from May 2000 to December 2002. The database allows to reconstruct the 
order book and thus there is no ambiguity in the data, especially for identifying 
the initiator of each trade. We present here the data of two stocks, namely 
Astrazeneca (symbol AZN) and Vodafone (symbol VOD). 

The second dataset contains a small sample of two NASDAQ stocks, namely 
Apple (AAPL) and Amazon (AMZN), for the period July- August 2009. The 
database contains trades and quotes and in this case we infer the initiator of the 
trades by using the Lee and Ready algorithm 17 . The algorithm consists in 

The equation found in [2] is slightly different because they define the temporary impact 
and the variance as proportional to the interval length r = T/N. 
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comparing the trade price with the prevailing quote and in identifying the trade 
as buyer (seller) initiated if the trade price is higher (lower) than the mid price 
of the prevailing quote. The case where trade price equals mid price remains 
undetermined, but it represents a very small fraction of the total (smaller than 
0.1%). Table [j] shows some summary information of the selected stocks. 





Number 


Avg. num. of 


Avg. price 


Avg. tick 


Avg. tick size - 


Symbol 


of days 


trans, in a day 


(£ or $) 


size (£ or $) 


price ratio (bp) 


AZN 


675 


731 


26.65 


0.01 


3.7 


VOD 


675 


1360 


2.95 


0.0025 


8.4 


AAPL 


42 


24094 


157.06 


0.01 


0.64 


AMZN 


42 


13890 


83.87 


0.01 


1.2 



Table 1: Summary information of the investigated stocks. 



4 Propagator model in trade time 

In order to solve the problem of the optimal execution, it is necessary to choose 
a model for the price dynamics. The most delicate aspect is to model in a 
careful way how price reacts to our trades, i.e. we need to have a good model 
of market impact. Here we propose to adapt a market impact model recently 
introduced [1] and termed propagator model. This model takes into account 
both the volume and the temporal dependence of the impact. In fact, it has 
been shown that price return reacts to a trade not only simultaneously to the 
trade, but also later in the future, and this effect is slowly decaying. In the next 
section we will review the propagator impact model. 

It is important to stress that the propagator model has been proposed and 
calibrated on transaction by transaction data. In the optimal execution scheme 
presented in the previous section, the length r of the time intervals corresponds 
to several minutes or the corresponding average number of trades. Therefore 
we need to adapt the model to describe impact in macroscopic time and test it 
on the data. This will be done in the next Section. 

The propagator model for market impact describes the price dynamics in 
transaction time, i.e. time t n , n € M increases of one unit after each transaction. 
The price p n is the log mid-price of the asset right before the transaction at time 
t n . Let v n be the signed volume of trade at time t n . 

The propagator model describes the price dynamics as 

n-l 

Pn = Po + X! fc ' v k) + Vk] ( 14 ) 

fc=0 

where Q is a general impact function which describes the impact at time t n of the 
trade of volume at time tk , r\k is an independent and identically distributed 
random term with mean and variance a 2 and po is the log mid price at a time 
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to far in the pastj^J In order to ensure causality G(J,k,Vk) is nonzero only for 
j > k. Moreover, Refs. |JJ showed from real data that the volume and time 
dependence of Q can beapproximately factoriscd, i.e. 

g{j,k,v k )^f(v k )-G {j,k). (15) 

Finally, we assume that Go is translationally invariant PQ, i.e. 

G (j,k) = G (j-k). (16) 

Since Q is non-zero only for j > k we have that Go(fc) is non-zero only for k > 0. 

Thus the model of price dynamics is a geometric random walk with market 
impact, i.e. 

j-i 

Pj =Po + $>* + f(v k )G (j - k)], (17) 

fc=0 

The difference between two consecutive log prices (i.e. the log-return) is 

rj =p j+1 - Pj =Go(l)f(v j )+J2[Go(k + l)-G (k)}f(v j . k ) + ri j . (18) 

fc>0 

The return rj can be divided into three components: the impact of the last 
transaction, which is proportional to f(yj), the effect of the decay of the impacts 
of all older transactions oc f{vj- k ), and a random shock r/j. 

The propagator Go of real transaction by transaction data is well fitted by 
the function [JJ ITS] 

Go(0= (zg + zV 2 ' (19) 

where Tq is a multiplicative constant, Zo is a correction for small lags and /3 the 
coefficient of the power law decay for Z ^> Z . 

The factorization of Q in a volume and a time component leads to the def- 
inition of the instantaneous impact function f(v). This function is the tick by 
tick expected return due to a trade initiated by a given signed volume v, i.e. 

f(v) = E[r\v], (20) 

Data indicates that f(v) is well approximated by an odd function, i.e. f(—v) ~ 
~f(v). 

Moreover, many studies [JJ, [TU], [5] have shown that / is a concave function 
of the volume v and is typically well fitted by a logarithmic function or a power 
law function with exponent smaller than one. Note however that this is true 
for individual trades returns. As returns are computed on longer time scales, 
the return conditional to volume imbalance becomes more and more linear (see 
Ref. [5]). 

2 In the original formulation this time was set at — oo. Here we consider a generic initial 
time (often set as the beginning of the trading program). 
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5 The propagator model under other time mea- 



sures 



The propagator model was originally introduced to describe the transaction 
by transaction dynamics of price. However, in order to use this model for 
the optimal execution problem as stated above, we need to have a model in 
aggregate time, i.e. where time advances by one unit after one interval r. First 
of all, notice that there are several possible definitions of aggregated time in 
financial data. The first, and most natural one, is obtained by considering r as 
a fixed real time interval, for example five minutes. This choice brings sometime 
problems related to the fact that financial processes are often not homogeneous 
in real time, because of periodicities, trading clustering, etc. For this reason one 
can consider also other time measures, such as aggregated trade (or transaction) 
time and aggregated volume time. In the former case, time advances by one unit 
any time a certain fixed number of trades have been executed. In the second 
measure, time advances by one unit when a given volume has been traded in 
the market. In this section, we will show how to adapt the propagator model 
to real time. In the Appendix we present the case of aggregated trade time. 



5.1 Calibration of the propagator model in real time 

We fix to as the opening time of the market and we define the series of times 
t n as separated by a constant time interval. For definiteness in the following 
we will consider 5 minute intervals, even if similar results apply to different 
durations. In this case to = 8:00, t x = 8:05, t 2 = 8:10, ... As above, p n is the log 
mid price right before time t n . We define the scries of aggregated volumes v n 
in terms of the volumes w* of the single transactions (tt stands for transaction 
by transaction), i.e. 

Vn= £ Vf. (21) 

[t„,t»+l] 



by this notation we mean that we sum the all signed volumes, i>'*, that are 
traded between time t n and i„+i. The quantity v n is the volume imbalance 
over the n-th time interval. This quantity may be very big in absolute value 
during time intervals of high market activity, like, for example, near market 
opening and closing. Therefore, we consider a related variable, the normalized 
volume imbalance w" or , which is obtained by dividing the volume imbalance v n 
in the n-th interval by the total absolute volume traded in the interval, 

„.nor "[tmtn+l] * (22) 



^[t„,*„+i] K I 

When there is no trade in the n-th interval, we define v™ or — 0. Note that 
u" or G [—1,1]- This definition has also been used in Ref. [10]. We can finally 
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state the propagator model in the new time setting: 

j'-i 

Pi=Po+J2l r )k + m° r )Go(j-k)}, n~IID (23) 
The difference between two consecutive log prices (log-return) is 



r, = p k+1 - Pk = J2 + H- ( 24 ) 

k=0 

where we defined g(k) = Go(fc + 1) — Go(fc), and Go(0) = 0. 

We now calibrate the propagator model on the four investigated stocks. We 
first estimate the impact function f(v nor ) of the normalized volume imbalance. 
It is defined as the expected return in a five minutes interval conditioned to a 
normalized volume imbalance v™ or in the same interval, i.e. 

f(v nor ) = E[r n \v™ r ). (25) 

In figure [I] we show the scatter plots of f(v nor ). As it is immediately evident, 
on this time scale and having used the normalization procedure, the impact 
function becomes very close to a linear function of the normalized volume. This 
will be a key ingredient to find the optimal execution, because the linearity of 
the impact allows an analytical solution of the optimization problem. For this 
reason we fit the impact function in real time with a linear function given by 

f(v nor ) =0v nor . (26) 

Table [2] reports the results of the fit for the four stocks. 

We can now estimate the impact propagator G in real time, by using the 



linear function f(v nor ) obtained above. We use eq. 24 as a linear regression 
model to find the best coefficients g(k) and from those obtain Go- From figure 
[2] we observe that for both datasets the propagator model seems to fit quite well 
the data. The fit is quite good also for the NASDAQ stocks, where the sample is 
small. This encourages to think that one does not need long historical datasets 
in order to fit the propagator model. 

We perform a best fit of the propagator with the functional form of Eq. 
19 The parameters of the best fits are reported in Table [2] Below we will 



use impact propagators with these parameters in order to calibrate the optimal 
trading schedules. 



5.2 Variance explained by the model 

We conclude this section by discussing more quantitatively how well the prop- 
agator model fits the data for different time definitions and time intervals. In 
particular, for each regression we compute the R 2 <G [0, 1] coefficient, which is 
the fraction of the total variance explained by the model. We summarize in 
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VOD - 5 minute aggreg. 
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AZN - 5 minute aggreg. 




5 minutes normalized ii 



5 minutes normalized imbalance 



AMZN - 5 minute aggreg. 




AAPL - 5 minute aggreg. 



5 minutes normalized imbalance 



5 minutes normalized ii 



Figure 1: Estimation of the normalized impact of Eq. 25 in real time. Parame- 
ters of the fits are given in Table [2]. 









Go parameters 




Symbol 


9 (bp) 


To 


lo 


P 


AZN 
VOD 
AAPL 
AMZN 


15.4 ±0.5 
26.0 ±0.8 
21.9 ±0.9 
26.9 ±0.8 


1.40 ±0.04 
1.07 ±0.01 
1.01 ±0.02 
1.05 ±0.03 


20 ± 1 

4± 1 
0.41 ±0.27 
0.70 ±0.35 


0.190 ±0.006 
0.075 ± 0.002 
0.23 ±0.007 
0.23 ±0.011 



Table 2: Parameters of the fits of the market impact propagator model in the 
real time setting, which we use in the impact model for the optimal execution. 
The parameters are those obtained from the fits shown in figure [T] and [2j 
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20 40 60 80 100 120 

5 minutes intervals 




5 minutes intervals 





Figure 2: Estimation of the impact propagator Go in real time five minute 
intervals. Each point is the result of binning many transactions with similar 
normalized imbalance and the y axis is the expected return. Parameters of the 
fits are given in Table [2j 
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Table [3] the R 2 coefficient for all the studied models. Apart from the cases 
discussed above, we include also the propagator fitted in trade time (as in the 
original paper of Bouchaud et al.), for aggregated transaction time with d = 8 
and d = 64 transactions, and for real time with 1 minute intervals. 







R 2 coefficients of the linear 


regression 




Symbol 


trade 


8 aggregated 


64 aggregated 


1 minute 


5 minutes 




time 


transactions 


transactions a 


ggregation 


aggregation 


AZN 


0.224 


0.245 


0.304 


0.192 


0.204 


VOD 


0.257 


0.329 


0.440 


0.209 


0.292 


AMZN 


0.172 


0.253 


0.229 


0.163 


0.153 


AAPL 


0.202 


0.302 


0.319 


0.223 


0.212 



Table 3: R 2 coefficients of the regression of the propagator model for different 
time definitions. 



We note that R 2 increases by increasing the length of the time interval (in 
number of transactions or in real time) . As an extreme case note that for VOD 
in the 64 aggregated transactions model, 44% of the total variance is explained 
by the model. Also the NASDAQ stocks are quite well fit by the model, even 
if when aggregation increases, the quality of the fit remains roughly constant, 
probably because of the smallness of the sample. 

6 Optimal execution model 

We now present the theoretical solution of the optimal execution in the presence 
of transient impact. For the case where bid ask spread costs are neglected we re- 
obtain the results presented recently in Ref. [5J [T^] . We define the logarithmic 
transaction costs c(v) , where pk is the log of the effective price P& 

N-l 

c(v) = ^2 Vk (P k ~ Po)- ( 27 ) 

fc=0 

We note that in the limit of short time horizons of the execution T, the 
effective price does not deviate much from the starting value Po. We can 
thus approximate: 

The logarithmic execution cost c(v) can therefore be thought as a fractional ex- 
ecution cost. More important, minimizing the execution cost C(v) is equivalent 
to minimizing the logarithmic execution cost c(v). 
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6.1 Baseline case: risk neutral investor without spread 
costs 

We are now ready to express the execution costs within the framework of the 
market impact propagator model described above. We start from equation |23[ 
which we report here (dropping the superscript nor of vJ? or ): 

n-l 

Pn = PO 

fe=0 

We recall that the effective log-midprice pk is the logarithm of the average 
mid-price at which we assume to trade the shares Vk between time tk and time 
tk+\. The log price p of the stock is a continuous time process, which we sample 
every 5 minutes to obtain the discrete price series p n . Therefore, it is reasonable 
to assume that the effective log midprice pk of the trade between time tk and 
tk+i is given by the average of the prices at the two times: 

Pk+Pk + l ,„„s 

Pk = g ' ( 30 > 

The equation that describes the dynamics of effective price is therefore 

n 

Pn = Po + ^[Vk + f{vk)G (n - k)] (31) 

fc=0 

where we defined the effective propagator Go as 

r fn\ G °W r m G "( 1 ) + G "( 2 ) r ^ G (2) + G (3) 

G o(UJ - — - — , G (1J - g ' G o(^J - ^ ' ••• ^ 6Z > 

With this impact model the fractional execution costs c(v) is given by 



N-l 

c(v) = V " 
n=0 

and its expected value is 



Tb 



.k=0 



(33) 



N-l 



E[c(v)] = <J2 f(vk)G (n - k)} (34) 

n=0 fe=0 

Let us assume for simplicity that that / is linear, even if it can be time 
dependent, i.e. 

f{vk)=O k v k . (35) 

where Ok are coefficients describing a possible deterministic change of impact 
during the trading day. As we have seen above, for five minute intervals the 
linearity of / is a reasonable approximation. 
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We now define the N x N triangular impact matrix as 

^ = {f 0ii - j) (36) 
which allows us to rewrite the expected logarithmic execution costs as 

E[c(v)} = v r Iv. (37) 
We want to find the optimal trading schedule, i.e. the one that minimizes 

minv T Zv, (38) 

V 

subject to the constraint 

Vk = l T v = X. (39) 

fc 

where 1 is a vector whose N elements are all 1. This a quadratic minimization 
problem, analogous to the well known Markowitz best portfolio problem (see 
[El HE]). We use a Lagrange multiplier z to enforce the constraint and we 
obtain the following unconstrained minimization problem 

min A(v, z) = min(v T Iv + z ■ (l T v - X)). (40) 

V V 

The optimal solution v* is 

v = -| r'l-j^r'i, (4i) 

and its expected execution costs: 

£[c(v*)] = v-Iv* = = I ^ ir (42) 

We note that the execution costs are quadratic in the total volume X. In- 
terestingly the optimal trading trajectory v* can be computed with a simple 
closed-form matrix expression. 

6.2 Optimal execution with spread costs 



One important extension to the model of last section is to consider the contribu- 
tion of bid ask spread costs. We assume that the bid-ask spread has a constant 
values — A — B, where A is the ask and B the bid. Recalling that P = A ^ B 
is the stock mid price, we define 5 as half of the fractional bid-ask spread, i.e. 

_ s/2 A-B 
S =^ = A^B- (43) 



3 The model can easily extended to the case of a deterministically varying (average) spread. 
This extension could for example capture the intraday variation of spread. 
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Then we express the effective log-price pk is 



Pn 



Po + / Im + f{vk)G (n ~ k)] + sign(vk)S, 



k=0 



(44) 



because we pay half the bid-ask spread on execution. 

The expected value of the transaction costs including the bid-ask spread 
contribution is therefore 



E[c(v)} 



v T Zv + <Sl T |v| 



(45) 



Minimizing Eq. 45 is more complicated than minimizing Eq. 37 because of the 
term with the absolute value of v. We can not use anymore the machinery of 
Lagrange multipliers and matrix inverse and we will use numerical minimization 
methods. 



6.3 Optimal execution with risk aversion 



In this section we consider an extension of the optimal execution model taking 
into account risk-aversion. The idea that risk matters in determining the optimal 
execution schedule was first proposed by Almgren & Chriss in [4]. We express 
the variance Var[c(v)] of the fractional transaction costs c(v) of equation 33 
Recalling that the rj n are independent and identical distributed random variables 
with zero mean and variance a 2 , we have 



N-l 



k-1 



Var[c(v)] 



(c(v)-£[c(v)]) 2 =£[(5> fc J>) 2 ] 



fc=0 j=o 
N-l N-l 



N-l N-l 

= £[(E^E^) 2 ]^ 2 E(E^) 2 - ( 46 ) 

fc=0 j=k k=0 j = k 

The variance is bilinear in the trading schedule v. We can therefore express 
it as 

Var[c(y)]=Y f VkjV k v j , (47) 

where V is 



(48) 



We now formulate the problem of minimizing execution costs with a risk- 
aversion term by following the construction of [4]. Specifically, we add the 
variance term multiplied by a coefficient A (which represents the risk aversion 
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parameter) to the transaction costs of equation 45 We obtain a total costs 
function c\(v) 

c A (v) = E[c(v t )] + XVar[c(v t )} + S(w) = v T [I + AV]v + 61 T \v\. (49) 

The minimum of the objective function in equation [49] can be found using nu- 
merical methods. 

Note that if spread costs are negligible we can solve analytically the opti- 
mization problem and find the solution that minimizes the objective function 



49 Setting S = 0, we define T = X + XV. By using again Lagrange multipliers, 



we have therefore the optimal trading schedule 

v* = z T-H = -^-^T-'l. (50) 

7 Empirically calibrated optimal execution strate- 
gies 

In this section we show a calibration to real data of the optimal execution 
schedules shown above. We will consider real time aggregation over five minute 
intervals. 

In section |5.1| we showed how to fit the propagator model when using real 
time intervals. Here we point out that some care should be taken when defining 
the impact function to be used in the optimization. In fact, the impact function 

f«° r ) = E[r n K or ] = <° r - (51) 

gives the relation between price change and normalized volume imbalance. We 
assume that the volume traded in the market has on average sign 0, i.e. trades 
are equally likely buyer and seller initiated. We define the total volume Wk 
traded between time tk and time tk+i as 

W k = J2 K"l ( 52 ) 

[tk i*fc+l] 

where vf refer to the volumes of the individual transactions. This quantity 



equals the denominator of the equation 22 The volume we trade at step k of 
the trading schedule gives a normalized imbalance of v^ or = Vk/Wk (because 
all the other volume has on average sign 0). Therefore the impact function / of 
the volume Vk is 

f(vD = ^ F Vk^e k Vk. (53) 
Wk 

which defines the series Ok- In the following we will approximate the series Wk 
with a constant value, Wk = W . In principle one could relax this assumption 
and take into account the intraday pattern of volume. The parameter 9 was 
estimated empirically and we report its value in Table [2] 
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In the following we will consider an investor that wants to execute a buy 
program in a trading day. We assume that the investor wants to execute a 
given fraction of the daily volume, and for definiteness we set this fraction at 
1%. We define the participation rate x as the ./V-dimensional vector, whose 
component Xk is 

x k =— . (54) 

i.e. the ratio between the traded volume Vk and the total market volume traded 
in the same interval, Wk- Given the 1% volume target and the approximation 
Wk — W, we have 



N-l , N-l 

X = 



N - 1 1 N - 1 AT Txr 

— F^ = " (55) 



100 100 

fc=0 fc=0 



Hence the participation rate x satisfies 

N-l 



N 

loo 



(a*) = 1% (56) 



fc=0 

We use the functional form of Eq. [19] for the impact propagator Gq. 



7.1 Optimal execution without spread costs 

We start with the simplest specification of the optimization problem by ne- 
glecting bid-ask spread costs and risk aversion. We computed analytically the 
optimal solutions v* for the four stocks. In figure [3] we plot the participation 
rate x*. 

We see that the participation rate of figure [3] oscillates between positive 
and negative values, but its average value is indeed 1%. The optimal solution 
consists in alternately buying and selling shares. When an investor buys he 
drives the price up. Thus, if he sells right away he takes advantage of the higher 
price and moves it back downj^] 

However, the solutions of figure [3] are of little practical interest. In reality 
one does not trade at the mid-price of the stock, but buys at the ask price 
and sells at the bid (because we are considering execution with market orders). 
By buying, one therefore pays the mid-price plus half the bid-ask spread. By 
selling, one gets the mid-price minus half the bid-ask spread. Thus, alternately 
buying and selling bears substantial costs of bid-ask spread. We show in the 
next section how including the bid-ask spread costs changes dramatically the 
solution. 

4 We note that the total impact costs are in any case positive, so that one does not profit 
from market impact, as studied in Ref. |19| . 
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VOD - 1% of daily volume 




AZN - 1% of daily volume 



3:05 09:45 11:25 13:05 14:45 16:25 8:05 09:45 11:25 13:05 14:45 lfi:25 




Figure 3: Optimal execution strategies to buy 1% of daily volume, neglecting the 
contribution of bid-ask spread and without risk aversion. We plot the optimal 
participation rate over each 5 minutes interval. 
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7.2 Optimal execution with spread costs 

We now solve the problem of optimal execution by including the cost of bid-ask 
spread. We approximate the bid ask spread with a constant value 25 throughout 
the day, neglecting for simplicity the intraday pattern of spread. The expected 
value of the fractional transaction costs is 

E[c(v)] = v T Xv + <Sl T |v|, (57) 

where 1 is the vector whose elements are all 1. The impact matrix X is the same 
as in section |7.1| and we estimated the parameter 6 from the data: 

$azn = 5.27 bp, Svod = 10.12 bp, 5 AM zn = 1.47bp, Saapl = 0.526p 

. (58) 

As we noted in section [6T2| we cannot express in a closed form the solution v* 
that minimizes Eq. |45| We instead applied a numerical optimization method, 
using the algorithm L-BFGS-B (see [20]: IS]) from the optimization routines 
of SciPy [55] . The resulting optimal participation rates x* are shown in Figure 

H 

We see that including the bid-ask spread contribution in the execution costs 
has radically changed the form of the optimal solution. The participation rate 
x* (and thus the trading schedule) is positive over all the market day. In other 
words, if one has to pay a bid-ask spread cost on every transaction it is not 
convenient to sell and buy back shares. We can see that the optimal solutions 
are U-shaped, i.e participation rate is high near market opening and closing, 
and low during the central hours of the day. This happens because the market 
impact propagator Go decays in time. In fact, the high participation rate at the 
opening impacts much the price. Then the period of low activity lets the market 
absorb the impact and recover a more favorable price. The high participation 
rate at the end instead impacts the price in the future, outside the time horizon 
T of the execution. 

The regularization achieved by the bid ask term is similar to what happens in 
portfolio optimization. It is known that adding to Markowitz objective function 
a penalty proportional to the sum of the absolute values of the portfolio weights 
stabilizes the solution and corresponds to an exclusion of short positions (23l 
124] . Similarly here, the bid ask spread term (which is motivated by the type 
of considered executions) stabilizes the solutions of the optimal execution and 
excludes oscillating solutions. We have tested that by choosing a S parameter 
much smaller than the fractional spread, one recovers the U-shaped solution. 

We compare the optimal solution v* just obtained with the flat solution of 
Bertsimas and Lo [3]. It consists of constant Vk — X/N over each interval. 
In Table [4] we show a summary of the execution costs of the optimal solutions 
analyzed, in particular the fractional (per share) spread costs c sp (v) and the 
fractional impact costs Ci mp (v), defined as 

(Sl T |v| v T Iv 
c *p( v ) = hT v i » <w(v) = (59) 
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We show these costs for all four stocks and three different solutions, namely 
the flat solution Eq. |8]of [HE], the U-shaped solution and the oscillating solution 
of section 17.11 







Execution costs per share, 1% averaj 


;e daily volume 




Fractional impact costs (bp) 


Fractional spread costs (bp) 




flat 


U-shaped 


oscillating 


flat 


U-shaped 


oscillating 


Symbol 


solution solution 


solution 


solution 


solution 


solution 


AZN 


4.36 


4.29 


4.20 


5.27 


5.27 


73.52 


VOD 


9.82 


9.76 


9.67 


10.12 


10.12 


81.62 


AMZN 


4.09 


4.03 


4.02 


1.47 


1.47 


14.20 


AAPL 


3.17 


3.12 


3.12 


0.52 


0.52 


6.15 



Table 4: Comparison of the fractional cost of impact and fractional cost of 
spread (see Eq. 59 ) for three solutions of the optimization problem, namely the 
flat solution of Bcrtsimas and Lo [3] , the U-shaped solution and the oscillating 
solution. 



We can see that if we are concerned with the minimization of the impact 
costs alone, the oscillating solutions are indeed the best ones. However, when 
we consider the contribution of bid-ask spread, the oscillating solutions become 
much less attractive: by selling shares and buying them back one incurs in 
very high spread costs. Instead, the flat and U-shaped solutions have fractional 
spread costs per share equal to S (because all elements of v have the same sign). 
The comparison between the U-shaped solution and the flat one shows that the 
former has impact costs between 1 and 2% lower. 



7.3 Risk aversion and efficient frontier of optimal execu- 
tion 

We now consider optimal execution by including a risk aversion term in the 
objective function. The objective function is (see Eq. 49 1 

ct t(v) = v T [Z + AV]v + 61 T |v|. (60) 

To express the variance matrix V we need to estimate the variance a 2 of the 
random shocks rj n . We calibrate it empirically as: 

* 2 =(i2vi)/L, (61) 

\n=0 / 
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where L is the total number of intervals in the dataset. The resulting average 
variances a 2 for the four stocks ar^El 

o-azn = 350.81 (bp) 2 , a 2 VOD = 764.52 (bp) 2 , 

vamzn = 395.62 (bp) 2 , a\ APL = 195.95 (bp) 2 . 



AMZN, Lambda = 0.0 AMZN, Lambda = 0.8 

2.5. 1 1 1 1 1 1 ^| 6 1 1 1 1 ^- 




AMZN, Lambda = 1.6 AMZN, Lambda = 2.4 




Figure 5: Optimal execution strategy to buy 1% of daily volume of AMZN for 
four different values of the risk aversion parameter A. 

We apply the same numerical optimization routine we described in the last 
section, for the minimization of the objective function of Eq. |60| The result 
for the stock AMZN is shown in Figure [5] Specifically, we plot the optimal 
participation rate profile x* for different values of the risk aversion parameter 
A. As the coefficient A increases, the U-shaped solutions becomes more and 

5 We note that the variance a 2 has an intraday profile: there are intervals (e.g. the ones 
near the closing time) whose returns have statistically higher variance than the returns over 
other intervals. We will not model this effect. 
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more front loaded, so that more trading activity is concentrated at the market 
opening rather than later in the day. This happens because by trading sooner 
one is less exposed to the fluctuations of the price and therefore the risk is lower. 

As a comparison we consider the optimal solutions with risk-aversion pro- 
posed by Almgren and Chriss [I] (see Eq. 12). 

In the case A = 0, Almgren and Chriss' solution is flat as in the Bertismas Lo 
model. We already know that the A = optimal solutions of propagator model 
outperforms the flat solution by 1-2% in impact costs. In order to compare the 
two models when A > 0, we plot the efficient frontier of the optimal execution 
according to the two models. Specifically, for each value of A we compute the 
expected cost and its variance and we plot them on a cartesian plane. The result 
is shown in figure [6] The dots correspond to Almgren and Chriss' solutions, 
the stars to the solutions resulting from our minimization procedure. We see 
that our optimal solutions outperform the v^ c solutions on the whole range of 
parameters of risk aversion A. For any given level of the variance of the cost, our 
optimal solution has lower expected impact costs than the corresponding iy& 
of an amount that is roughly constant (at about 1-2%). It is important to stress 
however that this is an in sample exercise and therefore the results might be 
quite different if one considers out of sample real executions. This comparison 
serves mainly to give an idea of the difference in the results obtained with the 
two execution strategies. 



8 Conclusions 

In conclusion we have presented an in depth analysis of the way in which the 
transient impact model introduced and discussed in Refs. [TJ [5] to describe the 
relation between order flow and price at the tick by tick level can be extended 
to describe impact on a longer time scale. We have considered (i) real time 
and (ii) aggregated trade time, where time advances by one unit any time a 
predetermined number of trades is executed in the market (see Appendix). We 
have shown how to calibrate the model and our empirical analysis has shown that 
the model fits quite well real data in both settings. This good agreement between 
model and data opens up the possibility of using the transient impact model for 
the optimal execution of large orders. By focusing on the real time setting, we 
have shown how the the model can be used to calibrate optimal execution in 
different cases. Specifically we have considered the role of risk aversion (following 
the approach of Ref. [1]), and of bid ask spread costs. Interestingly, when spread 
costs are neglected the solution turns out to be composed of alternating buy and 
sell trading intervals, irrespectively of the position one wants to take with the 
trade. The introduction of bid ask spread costs regularizes the solution, i.e. one 
always buys for buy trades. Numerical in sample calculations show that the 
proposed execution algorithm outperform existing algorithms. 



25 



VOD - 1% of daily volume 



Almgren & Chriss' solutions 
Results of the minimization 



AZN - 1% of daily volume 



Almgren & Chriss' solutions 
Results of the minimization 
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4 2^I00 3000 4000 5000 6000 7000 8000 9000 1000011000 
Variance of execution costs (bp)" 



AMZN - 1% of daily volume 



Almgren & Chriss' solutions 
Results of the minimization 



AAPL - 1% of daily volume 



Almgren & Chriss' solutions 
Results of the minimization 



2000 3000 4000 5000 6000 7000 8000 9000 10000 
Variance of execution costs (bp) 2 



2000 2500 3000 3500 4000 
Variance of execution costs (bp)* 



Figure 6: Optimal frontier to buy 1% of daily volume, for the four considered 
stocks. Stars correspond to the frontier obtained by using the propagator impact 
model, while dots refer to the frontier for the Almgren and Chriss model. On 
the x axis we have the variance of the execution costs, and on the y axis the 
expected value of the execution costs (per share). 
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9 Appendix: propagator model in aggregated 
trade time 



In transaction by transaction time we increase time by one unit after each 
transaction. In aggregated trade time wc aggregate d € IN transactions, so that 
time is increased by one unit after d transactions on the market. 

Let us denote £g M the time of the first transaction of the period we con- 
sider ("ati" stands for aggregated trade time). Then, every d transactions we 
increase time by one unit and we have that i^" is the time of the (n ■ d)-th 
transaction. The stock price p^" is defined as the log-midprice right before the 
(n ■ d)-th transaction. The volume v^ u is the algebraic sum of volumes of the d 
transactions between time i^" and time f^l_ x . We can express the volume v% a 
in term of single trades volumes vf (tt stands for transaction by transaction) as 



Clearly for d = 1 we obtain again the transaction by transaction time. 
9.1 Calibration of the propagator model 

We start by estimating the impact function f(v att ). This is the expected value 
of the log mid price change in an interval of d trades conditional to a given 
volume imbalance v att . This quantity has been studied in many papers on 
market impact (see Ref. [5] and references therein). 

We show in figure the plots of f(v att ) in aggregated trade time when we 
aggregate d = 8 trades for the two sets. We note that for the long LSE dataset 
the function is close to be linear for small volumes v att , and reaches a constant 
value for high values of v att . For the NASDAQ stocks the noise level is very high 
but a clear saturation of the impact for large volumes emerges. We therefore fit 
the impact with the arctangent function, as 



We can now estimate the impact propagator Go in aggregated trade time, 
mimicking the procedure used in real time. The results are shown in figureJH 
We also fitted a functional form for the impact propagator Go as in eq. [19 
We note that the propagator model seems to fit quite well the impact decay in 
aggregated trade time, even if we point out that there is no theoretical reason of 
why the propagator Go should preserve its functional shape under aggregation. 
To the best of our knowledge the theoretical problem of the description of the 
propagator model under aggregation of transaction is still open. 



d-(n+l)-l 




(62) 



f(v att ) = 6»arctan(p v att ). 



(63) 
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;ure 7: Impact of 8 aggregated transactions for the four considered stocks. 
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Figure 8: Estimation of the impact propagator Go in aggregated time of 8 
transactions for the four considered stocks. 
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Disclaimer 



This material (Material) is provided for your information only and does not constitute: (i) 
research or a product of the J. P. Morgan research department, (ii) an offer to sell, a solicitation 
of an offer to buy, or a recommendation for any investment product or strategy, or (iii) any 
investment, legal or tax advice. You are solely responsible for deciding whether any invest- 
ment product or strategy is appropriate for you based upon your investment goals, financial 
situation and tolerance for risk. J. P. Morgan disclaims all representations and warranties in 
the information contained in this Material, whether express or implied, including, without 
limitation, any warranty of satisfactory quality, completeness, accuracy, fitness for a partic- 
ular purpose or non-infringement. The information contained herein is as of the date and 
time referenced in the Material and J. P. Morgan does not undertake any obligation to update 
such information. All content, data, statements and other information are not warranted as to 
completeness or accuracy and are subject to change without notice. J. P. Morgan disclaims any 
responsibility or liability, whether in contract, tort (including, without limitation, negligence), 
equity or otherwise, for the quality, accuracy or completeness of the information contained in 
this Material, and for any reliance on, or uses to which, this Material, is put, and you are 
solely responsible for any use to which you put such information. Without limiting any of 
the foregoing, to the fullest extent permitted by applicable law, in no event shall J. P. Morgan 
have any liability for any special, punitive, indirect, or consequential damages (including lost 
profits or lost opportunity), in connection with the information contained in this Material, 
even if notified of the possibility of such damages. Any comments or statements made herein 
do not necessarily reflect those of J. P. Morgan, its subsidiaries or its affiliates. Any unau- 
thorized use, dissemination, distribution or copying of this Material, in whole or in part, is 
strictly prohibited. 

© 2012 JP Morgan Chase & Co. All rights reserved. J. P. Morgan is the global brand 
name for J. P. Morgan Chase & Co. and its subsidiaries and affiliates worldwide. J. P. Mor- 
gan Cazenove is a marketing name for the U.K. investment banking businesses and EMEA 
cash equities and equity research businesses of JPMorgan Chase & Co. and its subsidiaries, 
conducted primarily through J. P. Morgan Securities Ltd. (JPMSL). Execution services are 
offered through J. P. Morgan Securities LLC ( member of FINRA, NYSE and SIPC), JPMSL 
( member of the London Stock Exchange and is authorized and regulated by the Finan- 
cial Services Authority. Registered in England &: Wales No. 2711006. Registered Office 
125 London Wall, London EC2Y 5AJ), J. P. Morgan Securities (Asia Pacific) Limited (CE 
number AAJ321) ( regulated by the Hong Kong Monetary Authority and the Securities and 
Futures Commission in Hong Kong), and other investment banking affiliates and subsidiaries 
of JP Morgan in other jurisdictions worldwide and registered with local authorities as appro- 
priate. Please consult http://www.jpmorgan.com/pages/jpmorgan/investbk/global for more 
information. Prime brokerage services, including clearing and custody services, are offered by 
J. P. Morgan Clearing Corp. (member of FINRA, NYSE and SIPC). Product names, company 
names, and logos mentioned herein are trademarks or registered trademarks of their respective 
owners. 
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